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ABSTRACT 

We examine solitons in theories with heavy fermions. These "quantum" soli- 
tons differ dramatically from semi-classical (perturbative) solitons because fermion 
loop effects are important when the Yukawa coupling is strong. We focus on kinks 
in a (1 + l)-dimensional <$> theory coupled to fermions; a large- N expansion is 
employed to treat the Yukawa coupling g nonperturbatively. A local expression for 
the fermion vacuum energy is derived using the WKB approximation for the Dirac 
eigenvalues. We find that fermion loop corrections increase the energy of the kink 
and (for large g) decrease its size. For large g, the energy of the quantum kink 
is proportional to g, and its size scales as 1/g, unlike the classical kink; we argue 
that these features are generic to quantum solitons in theories with strong Yukawa 
couplings. We also discuss the possible instability of fermions to solitons. 
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1. Introduction 



Topological solitons, despite their inherently nonperturbative character, are 
typically studied semi-classically, that is, in a perturbative expansion in the cou- 
pling constants [1]. The first term in this expansion, the classical soliton, is 
the solution to a nonlinear classical field equation. This solution is nonpertur- 
bative because its energy diverges as the coupling constants — which parametrize 
the nonlinearity — vanish. Perturbative corrections to the soliton are important: 
they split the degeneracies of the classical solution resulting from Poincare and in- 
ternal symmetries, and project the solitons onto eigenstates of momentum, angular 
momentum, and charge. If the coupling constants are small, however, corrections 
to the shape and energy of the soliton are small, and the classical description of 
the soliton is essentially accurate. 

If the couplings are large, on the other hand, there is no reason to expect the 
quantum soliton states to resemble the classical solitons, at least quantitatively. In 
general, the strong coupling behavior of solitons in a quantum field theory is not 
well known. One notable exception is the sine-Gordon kink in (1 + 1) dimensions; 
because of the equivalence of the sine-Gordon theory to the massive Thirring model 
[2], the sine-Gordon kink at strong coupling becomes a weakly-coupled fermion in 
the Thirring model, which is well described by perturbation theory. 

In this paper we study strongly-coupled solitons more generally, when such a 
fortuitous equivalence does not arise. We focus in particular on solitons in the- 
ories with large Yukawa couplings. One motivation for doing so is the following. 
Fermions can acquire mass through a Yukawa coupling to a scalar field with non- 
vanishing vacuum expectation value. Solitons in such theories often carry (possi- 
bly fractional) fermion number. It has recurrently been suggested that when the 
Yukawa coupling is large such a soliton may have less energy than a fermion in a 
constant scalar field background; consequently, fermions may be unstable to the 
formation of solitons [3-11]. To determine whether this is so, however, one must 
know the form and energy of solitons in a strongly-coupled theory, which may differ 
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appreciably from classical solitons. Indeed, we expect fermion loop corrections to 
significantly affect the solitons when the Yukawa coupling is large. 

One means of studying a strongly-coupled Yukawa theory is through a large- 
TV expansion [5,10-13]. To leading order in 1/N, the theory can be solved for 
arbitrary values of the Yukawa coupling. This expansion captures some of the 
strong-coupling behavior of the theory, which one hopes is representative even 
when N is not large. To carry out this expansion, we introduce N fermion flavors 
and choose the iV-dependence of the couplings so that the theory has a sensible 
iV — > oo limit, with only fermion loops contributing to Green functions to leading 
order in 1/N. The total contribution of the fermion loops can be summed in closed 
form to give the exact large- N effective action 

S eff [0] =S[<j>] -iN\ogdet(ilp), (1.1) 

where S [0] is the classical scalar field action and Jp is the Dirac operator in the 
presence of the field 0. 

Solitons in this large- N theory are c-number configurations of the scalar fields; 
scalar field fluctuations are suppressed because scalar loops do not contribute to 
the effective action to leading order in 1/N. The shape of the large- N soliton 
differs from the classical soliton, however, since it extremizes not the classical 
action but the effective action (1.1). The fermion loop contribution significantly 
alters the form of the soliton when the Yukawa coupling is large. In this regime, 
where quantum effects are so important, the large- N soliton is truly a "quantum 
soliton." 

To determine the form of the quantum soliton, we need to know — iN log det (iJp) 
explicitly for an arbitrary scalar field configuration. One generally resorts to 
some local approximation, such as the gradient expansion [5, 14], accurate for 
slowly-varying configurations. The gradient expansion, however, breaks down for 
topological solitons in the theories that we are considering. Another approach 
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to computing the fermion loop contribution relies on the fact that for static soli- 
tons (iN/T) logdet (ilp) is just the energy of the "Dirac sea," the sum of negative 
eigenvalues of the Dirac equation in the soliton background [15]. Unfortunately, the 
Dirac eigenvalues must be numerically computed [16] for each separate background 
considered, rendering this approach inconvenient for a variational problem. 

In this paper, we propose a hybrid of the gradient expansion and eigenvalue 
sum methods. Following an idea of Wasson and Koonin [17], we use the WKB 
approximation to estimate the Dirac sea eigenvalues for an arbitrary static scalar 
field background. We then sum these to obtain a local expression for the fermion 
vacuum energy. Unlike the gradient expansion, this expression is finite for topo- 
logically nontrivial configurations. Using this WKB approximation, we extremize 
the effective action to find the form of the quantum soliton in the large- N theory. 

We illustrate this method on a well-known example, the kink of the (1 + 1)- 
dimensional 4 theory coupled to fermions. The classical kink is reviewed in sect. 2. 
In sect. 3, we derive the WKB approximation for the large-iV effective action in 
this theory. This result is used in sect. 4 to find the form of the quantum kink, 
which is contrasted to the classical kink. The question of fermion stability is also 
discussed. In sect. 5, we present our conclusions and discuss the features of the 
model that we expect are generic to strongly-coupled solitons. 

2. Classical Kinks 

We begin by recalling the form and quantum numbers of the classical kink 
[1]. The (l + l)-dimensional 4 theory coupled to iV flavors of fermion has the 
Lagrangian 

The iV-dependence of the parameters has been chosen so that this theory has a 
sensible N — > oo limit. If we rewrite as \/~N(p, the parameter N becomes an 
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overall scale, 



\ [d^f - \\ (y. 2 - v 2 ) 2 ] M - 9<P) V- (2-2) 
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In the vacuum state \<p\ = v, the scalar field has mass \/2\v and the fermion 
field mass gv. In two dimensions, v is dimensionless, the scalar self-coupling A has 
dimension 2, and the Yukawa coupling g dimension 1. It is convenient to substitute 
for A and g the parameters 



The parameter x c \ is proportional to the scalar field Compton wavelength (and, 
as we will see, the size of the classical kink), and will serve as the overall scale of 
length and energy in the theory. There are two dimensionless parameters, v and 
y, the latter being proportional to the ratio of fermion and scalar masses. 

The Lagrangian (2.2) gives rise to the field equations 

1 N - 

d 2 v + V 3 _ Xv 2 v = _ g _J2^ (2.4) 

i=l 

(i^-g(p)^ = 0. (2.5) 

The topologically nontrivial solutions of these equations give a "classical" descrip- 
tion of the soliton states in the Hilbert space, which is accurate when the quantum 
corrections are small. If we neglect the fermion source term, the scalar field equa- 
tion (2.4) has the well-known static kink solution 

VciO^) = v tanh ( — j , (2.6) 

which is the lowest energy state with topological charge [<p(oo) — ip{— oo)]/2v = 1. 
There is also an anti-kink solution which interpolates from v to — v, with topological 
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charge —1. The Dirac equation (2.5) in the kink background (2.6) has a self- 
conjugate zero mode solution 

M*) = ( lsecH T« )] "), 7 1 = . <ri ' (27) 

V / 7 — 2°3- 

The state with the zero mode occupied has the same energy as that with the zero 
mode unoccupied. Since there is a zero mode for each flavor i, the kink is 2^-fold 
degenerate. If n of the zero modes are occupied, the kink has fermion number 
n — ijV, which ranges from — \N to [18]. The anti-kink too has degeneracy 
2^, and fermion number ranging from —\N to \N . Although the fermion zero 
modes increase the degeneracy of the kink, their contribution to the source term in 
the scalar field equation (2.4) vanishes, so the kink (2.6) remains a solution even 
in the presence of fermions. The energy of the classical kink 

EM=^NyKv* = lj(£) (2.8) 

has no dependence on the Yukawa coupling g, even though the kink carries fermion 
number due to the zero modes. 

This classical picture leads to the fascinating possibility that, even if fermion 
number is conserved, "ordinary" fermions may be unstable to the formation of 
solitons carrying fermion number. A configuration consisting of a widely-separated 
kink and anti-kink, each carrying fermion number ^N, has zero topological charge, 
fermion number N, and energy |t> 2 (N/x c \). On the other hand, a set of iV widely- 
separated fermions in the vacuum background <p — v, a state which has the same 
quantum numbers, has energy Ngv = y(N/x c \). Thus, when y > |t> 2 , it is energet- 
ically favorable for a state of N fermions to coalesce onto a spontaneously created 
kink/anti-kink pair. Each kink acts as a kind of bound state of fermions. Even 
more surprising, in a theory with one flavor of fermion (N = 1), a single fermion 
could split into a kink/anti-kink pair, each with fermion number ^. 
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This putative instability occurs only when y is large, however, where the quan- 
tum corrections from fermion loops are important and the semi-classical approxi- 
mation breaks down. To determine whether fermions are truly unstable, one must 
compare their energy not with that of a classical kink, but of a "quantum kink," 
which includes the effects of quantum corrections. The quantum kink extremizes 
not the action but rather the effective action. In the next section we will derive a 
local expression for the effective action suitable for finding the quantum kink. 



Quantum solitons are field configurations that extremize the effective action, 
which includes quantum corrections. To find the form of quantum solitons, one 
needs an explicit local expression for the effective action. The familiar gradient 
expansion, however, diverges for topologically nontrivial configurations in (1 + 1)- 
dimensional 4 theory. In this section, we derive an alternative local approximation 
for the effective action that is finite for kinks. 

Since we are interested in the properties of solitons for large Yukawa coupling 
g, the effective action must be calculated nonperturbatively in g. This can be done 
by taking the number N of fermion flavors to be large, holding A, v, and g fixed, 
and calculating to leading order in 1/N. Scalar field fluctuations are subleading in 
1/iV, so only fermion loops contribute to the large- N effective action 



to tame the divergent contributions of the fermion determinant to the one- and 



3. Effective Action for Kinks 




iNlogdet (i$-g<p) (3.1) 



We have added the counterterm 



6C(<p) = AN(<p 2 




(3.2) 



7 



two-point functions, and the overall constant iN log det (i<^ — gv) to ensure that 
S e ff[if — v] — 0. The coefficient A is fixed by requiring the one-point function to 
vanish at ip — v, 



= 2Av + ig 



d 2 p 



tr 



(3.3) 



(2tt) 2 \$ - gvj ' 

so that v remains the minimum of the effective potential. With a cutoff A on the 
spatial momentum pi, eq. (3.3) gives 

A 



o 



p\ + g 2 v 2 



(3.4) 



This counterterm also renders finite the two-point function 

,2 



rlS(p) 



p=0 



= -2Aw 



2 r 



7T 



(3.5) 



where a — tp — v. Fermion loop contributions to all other Green functions are 
finite. 

We must write the effective action (3.1) in a more tractable form if we are to 
find the quantum kink explicitly. The gradient expansion [5, 14] 



(3.6) 



is a useful approximation for slowly-varying fields. The first term in this expansion 
is minus the effective potential 



V eS (<f) A 



N 4 
The term with two derivatives is 



( ^)VfVlng)-f!(,W). (3.7) 



N 



1 + 



127ry^ 



(3.8) 



At this point, we discover that the gradient expansion fails for topological soli- 
tons in this theory; any configuration (p(x) with unit topological charge must pass 
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(2) 

through ip — somewhere, at which point C e g((p), as well as higher order terms, 
diverges. This failure is quite general. For the gradient expansion to converge, field 
gradients must be small relative to gtp, the "local fermion mass." Since the latter 
vanishes at the core of solitons with fermion zero modes, the gradient expansion 
necessarily breaks down there, no matter how slowly varying the field. 

An alternative approach for a static scalar field background such as the kink 
is to express the effective action in terms of Dirac equation eigenvalues [15,4]. For 
time-independent (p(x), the effective action equals — E e g [ip] T, where T = f dt and 
E c ff [ip] is the energy of the configuration, 



E eS [<p] = E cl [ip] + Q 



(3.9) 



a sum of the classical energy 
E c i [<p] = — 



oo 








/d, 


)( 






-oo 





( &V\ 2 , 1 / 2 2V 



X 



(3.10) 



and the quantum correction, the fermion vacuum energy, 

Q[<p] 



iN iN 

— log det (^ -gp)-— log det (^ - gv) + 5E [ip] . 



(3.11) 



The first term in eq. (3.11) can be interpreted as the energy of the Dirac sea in the 
background (p(x). 

To write eq. (3.11) more explicitly, we observe that the Dirac equation (2.5) im- 
plies that the spinor components ip 1 = (Jpt^j obey the Schrodinger-type equations 

d 2 



2 

± 



-y V ± {z)-y +x cl e 
in a static background <p(z), where 



± 







(3.12) 



V±(z) 



-I T 



1 dp> 



(3.13) 



y* J yv dz 

We restrict ip(z) to configurations of unit topological charge that obey <p(— z) = 
—ip(z); the Schrodinger potentials V a (z) are then even, and the solutions tpa(z) 
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can be taken to be parity eigenstates. (Here o = ± labels the upper and lower 
spinor components.) Since |0(±oo)| = v, the potentials V a (z) vanish at ±00, 
so eq. (3.12) has a continuous spectrum of states classified by their asymptotic 
momentum, k = \J x 2 ^ 2 — y 2 , and their parity. The asymptotic forms of the 
continuum wavef unctions 



*±oc 



(3.14) 

Vv,odd(M) — * sin(kz ± 5<5 ff dd(fc))> 

serve to define the phase shifts 5 ajeven (k) [5<r,odd(&)] f° r the even (odd) parity 
states. If we put the system into a box, \z\ < \h, with periodic boundary con- 
ditions, eq. (3.14) implies that the allowed momenta satisfy k an L + 5 a (k) = 2nn. 
Eq. (3.12) may also have a series of discrete bound states with eigenvalues e 2 ai < 
y 2 /x 2 v Because the configuration has topological charge 1, the upper spinor com- 
ponent is guaranteed [18] to have a zero mode, e+ = 0. Thus, any configuration 
with unit topological charge can carry fermion quantum numbers. 

The difference of fermion loop contributions can be written as the shift of the 
Dirac sea energy [15] 

^ log det (i? -g<p)~ det - gv) = - \N ^ " e S) > ( 3 - 15 ) 

<y A 

where e G \ denotes the positive root of e 2 x , and are the Dirac eigenvalues in the 
constant configuration <f(x) = v. Eq. (3.15) may be separated into the sum over 
discrete eigenvalues 

£ di sc [<p] = -^EE ( e - - — ) ( 3 - 16 ) 

a i V X c\J 

and the sum over continuum eigenvalues 



£cont M = - JiV^E E \< k ™) - < k °h] » = ' /k2+y2 . (3.17) 

o n>0 parity 

Using k an L + 8„{k) = k^L = 2irn, and letting L — > 00, we can write eq. (3.17) as 
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[15] 



A A 
a CT 

where 5 = | (5 even + 5 odd ). The integral over k diverges as the momentum cutoff 
A is removed, but this divergence is cancelled by the counterterm energy 

oo A oo 

w--hw-£(3fwT?h($- 1 )- (319) 

— oo — oo 

The sum of eqs. (3.16), (3.18), and (3.19), 

Q M = E disc [<p] + E cont [<p] + 5E [<p] (3.20) 

is precisely the fermion vacuum energy (3.11). 

The expression (3.20) for the fermion vacuum energy is much more explicit 
than eq. (3.1), and can even be computed analytically for certain scalar field con- 
figurations [19]. For an arbitrary background, however, e a i and 5 a (k) must be 
computed numerically [16]. Wasson and Koonin [17] showed how to speed up the 
convergence of these "brute force" numerical calculations by employing the WKB 
approximation for the high momentum phase shifts, but the discrete eigenvalues 
and low momentum phase shifts must still be computed numerically for each sepa- 
rate field configuration. Thus, eq. (3.20) is still not very convenient for extremizing 
the effective action* 

Taking our cue from ref. [17], we adopt the WKB approximation for all the 
Dirac eigenvalues, both continuous and discrete, and use them in eq. (3.20) to 
obtain a local expression for the energy of an arbitrary scalar field configuration. 

* Campbell and Liao [4] were able to extrcmizc (3.9) using powerful inverse scattering 
methods, but only for the special case y = 1. 
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The resulting expression will be accurate for field configurations slowly varying on 
the scale of the fermion Compton wavelength, but unlike the gradient expansion, 
does not diverge for solitons. We will then use this approximate expression to find 
the form of the quantum kink in sect. 4. 

In the WKB approximation, the continuum eigenfunctions of eq. (3.12) are 

<e K ven(^ *) = COS ( £ k a (z')dz') , 

V 1 , (3-21) 

€!odd(^ z) = -== sin (f* k a {z')dz>) , k a {z) = Jk* - ylV a (z), 

vK(z) v 

whence the phase shift defined through eq. (3.14) is given by 

oo 

C B W= I dz[k„(z)-k], (3.22) 

— oo 

independent of parity. (We assume V a (z) < everywhere; this will be true if (fi(z) 
does not vary too rapidly.) Using the WKB phase shifts (3.22) in the integral 
(3.18) and adding the counterterm energy (3.19), we find 

oo 

*5SM + / dz[(i-^)-(^ + ^) 

— oo 

+ (1 + V+) log (l + V / -^) + (1 + V-) log (l + y^V- 

(3.23) 

We also need to approximate the sum over discrete eigenvalues (3.16). In the 
WKB approximation, the Schrodinger equation (3.12) has discrete eigenvalues e 
whenever w a (e), defined by 

oo 

Me) = I J dz k a (z)e(k 2 a (z)), k a (z) = ^x 2 cl e 2 -y 2 -y 2 V a (z), (3.24) 

— oo 

equals half an odd integer, wGZ + j. The number of discrete eigenstates is given 
by the integer closest to w a (y/x c \). We define e a (w) by inverting eq. (3.24) and 
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setting e a {w) = for < w < w a {0). The sum over discrete eigenvalues (3.16) in 
the WKB approximation is then written 







(2) 

and E^ c [ip\ is the remainder. The integral (3.27) may be rewritten 



^disc 



= - hJ2 de w ^ 



o 

y/zci 



-oo CT 



47T (x c i) ' 



(3.25) 



We separate this into two terms 

£dT M = E^ c [<p] + E%i [<p] , (3.26) 



where E^ c [ip] is the integral approximation of the sum (3.25) 



w<,{y/x c i) 

e£L if] = - h E / dt » (^m - ^) • < 3 - 27 ) 



/d,[(^+v^) + (l + V' + )>og(^g^ 

— OO ^ 

+ a + v.,v,(^ 

(3.28) 

Adding the contributions from the continuum (3.23) and discrete (3.26) states, we 
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obtain 



Q WKB [<p] 



y^/N_ 

47T V^cl 



dz 



1 - 



+ 



1 Up 



(p 1 dip . , 

^ + -/ log 

f z yv dz 



1 dy? 
2 + — TT 



1 dy?\ 
yi> dz j 



1 d<^ 
yv dz 



E 



(2) 
disc 



(3.29) 



for the fermion vacuum energy in the WKB approximation. 



When ip(z) is slowly varying on the scale of the fermion Compton wavelength, 
the number of discrete states w a {y/x c \) is large, the sum (3.25) is well approximated 
by the integral (3.27), and E^? [ip] is much smaller than [ip] . If we therefore 

(2) 

neglect E\^ M > ec l- (3-29) provides a completely explicit local expression for the 
energy of a static configuration 



E 



off 



+ 
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oo 












-oo 


i. 


{-) 


4tt 


\Xc\J 



1 AV\ , 1 ( 2 



2\2 
u ) 



d^ 



1 - 



+ 



1 /y? 



1 d^ 

yv dz 



log 



+ 



1 (ip 



1 d</A 
yv dz j 



1 dtp 
yv dz 



1 dip 

yv dz 
(3.30) 

For ip(z) constant, eq. (3.30) reduces to the effective potential (3.7). When ip(z) is 
not constant, eq. (3.30) yields a correction to the effective potential which, unlike 
the gradient expansion, does not diverge for configurations going through ip — 0. 

We conclude this section by comparing the WKB approximation of the fermion 
vacuum energy of the classical kink, <p c i(z) = v tanh(z), with the known exact re- 
sult. The WKB approximation should be accurate for y ^> 1, when <p c \(z) is 
slowly-varying relative to the fermion Compton wavelength. The Dirac equation 
can be solved analytically in the classical kink background. Using the result- 
ing eigenvalues, Chang and Yan [19] computed the exact fermion loop correction 
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(3.20) to the energy of the classical kink 



Q bd] 



-) 



A(y). 



(3.31) 



The function A(y) is given by a complicated integral, but for integer y it simplifies 
to [19] 



y-i 



A(y) 



7T 



+ I —\/2yn — n 2 + — \/y 2 — n 2 arctan 



n=l 



- 1 



yez. 

(3.32) 



Using the Euler-Maclaurin formula, we obtain the large y behavior of eq. (3.32) 
3 



A(y) 



2tt 8 



y 



(i \ °° (Ah ?\\\\ 

k — 1 



.0206 . . . 



(3.33) 

where f^fc are the Bernoulli numbers. The series defining (5 is asymptotic, so we 
only keep 4 or 5 terms in the sum. 

The WKB approximation is obtained by substituting (p c i(z) into eq. (3.29) and 
expanding for large y 



Q WKB M 



N 



3 

27 



(3.34) 



Using the WKB approximation for the discrete eigenvalues e a i together with the 
Euler-Maclaurin formula, we find the leading behavior of the remainder term 



42c M = §- 



2 ^ i3+ i7f-l)^ + ° ( 



(3.35) 



Thus the WKB approximation of the fermion vacuum energy is 
N 



Q WKB M 



3 
27 



(3.36) 



in agreement with eq. (3.33) to this accuracy. The y 2 term, of course, is just the 
contribution from the effective potential (3.7). The non-analytic subleading ^Jy 
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dependence cannot be seen in the gradient expansion, but is correctly given by the 
WKB approximation (3.29). 

Obviously, the coefficient obtained for the subleading ^fy dependence would be 
incorrect if we made the further approximation of dropping i^isc M > as was done 
in obtaining the local expression (3.30). Nonetheless, eq. (3.30) correctly gives the 
order of the subleading dependence. In general, it provides a useful estimate of the 
correction to the effective potential for a spatially-varying field. 



4. Quantum Kinks 

The quantum kink extremizes the effective action of the (1+1) -dimensional c/> 4 
theory. For small Yukawa coupling y, the effective action (3.1) differs only slightly 
from the classical action, so the quantum kink nearly coincides with the classical 
kink. When y is large, however, fermion loop corrections are important, and the 
quantum kink differs significantly from the classical kink. 

To find the explicit form of the quantum kink, we use the local approximation 
(3.30) for the energy E e $ [if] of a static scalar field configuration derived in sect. 3. 
The equation of motion for the quantum kink follows from extremizing eq. (3.30), 

l y? 2 \ dV 9 ,V A y 2 , | i <p /<V 2 

(f 4 If dip^ 2 



./V.) — ■> ■> \ , 

(4.1) 

Using the program COLSYS [20], we have solved this equation numerically for 
various values of the parameters subject to the boundary condition ip(±oo) = 
±v. The solutions obtained interpolate smoothly between — v and v. Indeed, 
their profiles are almost indistinguishable from the hyperbolic tangent shape of 
the classical kink (see fig. 1). The slope of the quantum kink differs from that of 
the classical kink, however, being much steeper for certain values of the parameters. 
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We can more easily see how the slope of the quantum kink depends on the pa- 
rameters of the theory by restricting ip(x) to the one-parameter family of functions 



x 



¥x (x) = v tanh — , 



where xo is the "size" of the ansatz. We write the energy of the ansatz 

XQ 



(4.2) 



Ees(zo) = E c \(zq) + Q(zq), Zq = , 



(4.3) 



where z$ is the ratio of the size of the ansatz to that of the classical kink. The 
classical contribution 



E c \{zq) 



N 


"2 2 






-V 






3 





(4.4) 



has a minimum at zq — 1, of course. The quantum contribution is obtained by 
substituting the ansatz (4.2) into the WKB approximation (3.29) and retaining 
the leading power of y 



Q WKB (zo) 



N 



3 

27 



7T 



y 2 z + O 



(4.5) 



The WKB approximation is accurate when the neglected terms are small, which 
requires zq ^> 1/y. That is, the size of the ansatz must be much larger than the 
fermion Compton wavelength (xq ^> l/gv). 

In the following discussion, we assume large Yukawa coupling y ^> 1. The size 
of the quantum kink is found by minimizing (4.3), 



1 + 



9 3tt\ y 2 
47 ~ 16 J v 2 



-1/2 



(4.6) 



and depends on the values of both dimensionless parameters y and v. When v ^> y, 
the kink size zq ~ 1 and the quantum kink reduces to the classical kink, because the 
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classical contribution to the energy is dominant in this regime. On the other hand, 

_i 

when v <ti y (but v > 1), the kink size zq ~ (Jp — ff ) 2 (v/y) « 2.8(v/y); the 

quantum kink is much smaller than the classical kink. The energy of the quantum 

i 

kink in this regime, E e g (- — | ) 2 vy(N/x c \) ~ .48 vy(N/x c \), is larger than the 
classical kink energy (2.8) due to the positive fermion vacuum energy (4.5). 

When f < 1, the WKB approximation (4.5) breaks down because the kink size 
is no longer much larger than the fermion Compton wavelength. By using the exact 
Dirac eigenvalues for the background (4.2) rather than the WKB eigenvalues, how- 
ever, we can calculate the fermion vacuum energy Q{zq) without approximation, 
just as for the classical kink. We find 

Q(*o) = (—) yU(yz ), U(t) = M (4.7) 

\x c ij t 

where A(y) is the fermion vacuum energy of the classical kink defined in sect. 3. 
The function U(t) is shown in fig. 2, and equals l/ir at its minimum t — 1. (That 
its minimum is at t — 1 can be seen from eq. (3.32) and from 

dA 1 t £=i / r~^~ ' 2 ft^n ft 2 ~ \ 

at 2 7r \ \ 2t — n n \ t + n V r I 

(4.8) 

obtained by a calculation similar to that in ref. [19].) When t; C 1 (and y ^> 1), the 
fermion vacuum contribution (4.7) dominates the energy, so the size of the quantum 
kink is determined by the minimum of Q(zq), that is, zq ~ 1/y. The quantum kink 
energy, E e q xs (y/n)(N/x c i), is much larger than that of the classical kink (2.8). 

The classical contribution to the energy E c \{zq) is minimized for zq — 1, when 
the ansatz size equals the scalar field Compton wavelength, xq = x c \. The quantum 
contribution Q(zq) is minimized for zq = 1/y, when the ansatz size equals the 
fermion Compton wavelength, x$ = x c \/y = 1/gv. The size of the quantum kink 
always lies somewhere between these two values. The three limits we considered 
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above 



v < 1 < y =^ z 



1 < v < y 



1 < y < v =>- Zq 



1. 



t 

l 




(4.9) 



correspond to the regime in which the classical energy is dominant (v ^> y), the 
regime in which the fermion vacuum energy is dominant (t; C 1), and the regime in 
which both contributions are important (1< !) < y). For v > y, the classical and 
quantum kink nearly coincide, while for v < y, the quantum kink is smaller and 
has greater energy than the classical kink. Note that due to the fermion vacuum 
contribution (4.7), the energy of the kink is bounded below by (y/ir){N/x c \) = 
Ngv/n, that is, 1/ir times the mass of iV fermions. 

We now turn to the question of fermion stability. In sect. 2, we saw that 
for sufficiently strong Yukawa coupling, y > |t> 2 , a state of N widely-separated 
fermions has greater energy than a kink/anti-kink pair, computed in the classical 
approximation, so one might expect a kink/anti-kink pair to appear spontaneously, 
with the fermions coalescing to occupy the zero modes. Since the zero modes do 
not increase the kink energy, the energy of the fermions on the kinks is independent 
of y in the classical approximation, and would be much less than the energy of the 
fermions in a constant scalar field background for y ^> v 2 . The kink binding energy 
could approach 100% for very large Yukawa coupling. 

Instead we have found that, for large y, quantum corrections significantly in- 
crease the energy of the kink. For v > 1 (but v y), a kink/anti-kink pair 
has energy ~ yv(N/x c \), greater than the energy of N fermions, so the fermions 
are stable. For v < 1, the energy of a kink/anti-kink pair may be less than 
y(N/x c \) = Ngv, in which case a state of N fermions may be unstable to the for- 
mation of a kink and anti-kink, each carrying fermion number ^N. Since the kink 
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energy is never less than (y /ir)(N/x c i), however, the energy of a widely-separated 
kink/anti-kink pair is not significantly less than that of the original fermions; the 
binding energy per fermion cannot exceed 1 — ^ ~ 36%. 

Up to this point, we have been chiefly concerned with large Yukawa coupling, 
y >> 1; we conclude this section by briefly considering y < 1. When y is not large, 
the WKB approximation is no longer useful, but we can use the exact solution 
(4.7) for the ansatz (4.2). The case y — 1 is interesting, because then zq = 1 
minimizes both the classical and quantum contributions to the energy; the classical 
kink is an extremum of the effective action restricted to the subspace of functions 
(4.2). One might suspect from this that the classical kink extremizes the effective 
action over the space of all functions. Campbell and Liao [4] proved this to be the 
case by using inverse scattering methods (which were tractable only when y — 1). 
Thus, the quantum kink exactly coincides with the classical kink (for all values of 
v) when y = 1? As we have seen, they differ when y ^ 1. 

The energy of the kink when y — 1 is 

4 2 1\ N 



E eS = i-v z + - — , y = l, (4.10) 

so a kink/anti-kink pair will have less energy than N widely-separated fermions 



when v < y | (l — ^) .37. It turns out that iV fermions are unstable to kink 
formation only if v < \J~^ ~ -28; the most energetically favorable configuration of 
N fermions for v > is a bag [4]. (See also ref. [9].) 

Finally, for small Yukawa coupling, y < 1, the fermion Compton wavelength is 
larger than the scalar field Compton wavelength, so quantum corrections tend to 
increase the size of the kink. When C ?/ < 1, the quantum contribution domi- 
nates the energy, and the kink has size zq ps 1/y and energy E e g pa (y/n)(N/x c \). 



* The stationary phase approximation of ref. [4] is equivalent to our large- N approximation, 
f Interestingly, the theory is supersymmetric precisely when y = 1 [4, 21]. 
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5. Conclusions 



We have examined the effects of quantum corrections on solitons in a (1 + 1)- 
dimensional 4 theory with a large Yukawa coupling y to fermions. To treat the 
Yukawa coupling nonperturbatively, we have solved the theory in the large- N limit, 
where N is the number of flavors. The solitons in this theory are kinks which carry 
fermion number ranging from —\N to ^N. In the classical approximation, the 
energy of the kink is independent of y, and its size is proportional to the scalar 
field Compton wavelength. We have found that fermion loop corrections increase 
the energy of the kink and (when y > 1) reduce its size. As a result of the fermion 
vacuum contribution, the kink energy is bounded below by {y /ir)(N / x c \) = Ngv/n, 
and its size can be as small as the fermion Compton wavelength. 

When y is large, a state of N fermions is expected on classical grounds to be 
unstable to the formation of a kink and anti-kink, each carrying fermion number 
7;iV. Quantum corrections eliminate this instability for v > 1 by increasing the 
kink/anti-kink energy. The instability persists for v < 1, but the difference in 
energy between the N fermions and the kink/anti-kink pair is only about 36% 
because the kink energy is proportional to the Yukawa coupling in the large y 
limit. 

In the large- N limit, scalar loops are suppressed. The energy of scalar field fluc- 
tuations is of order l/x c \, small compared to the classical kink energy ~ Nv 2 /x c \. 
What happens when N is not large? Will a single fermion decay into a kink/anti- 
kink pair when N — 1? Scalar field fluctuations are still relatively unimportant as 
long as v is large; 1/v 2 is the usual semi-classical expansion parameter. We found 
fermions to be stable in this regime. Scalar corrections become more important 
for small v, but on the other hand l/x c \ is still small relative to the quantum kink 
energy y/ixx c \ when y is large. It is difficult to say whether a single fermion is 
unstable when v < 1. 

Some aspects of the model discussed in this paper are peculiar to two dimen- 
sions. Presumably only in two dimensions can a fermion split into a pair of solitons, 
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each carrying fermion number ^- We expect other features of the quantum kink to 
be more universal, however. First, its energy acquires a linear dependence on the 
Yukawa coupling in the strong coupling limit through the fermion vacuum energy 
Second, for large Yukawa coupling, fermion loop corrections tend to reduce the size 
of the soliton in the direction of the fermion Compton wavelength. Both of these 
features apply not only to the (l + l)-dimensional solitons described in this paper, 
but also to (3 + l)-dimensional (large- N) nontopological solitons [10]. 

General arguments can be adduced to suggest that these are generic features 
of quantum solitons in any (large- N) strongly-coupled Yukawa theory in (3+1) di- 
mensions. The fermion loop contribution to the effective action is — iN log det {ilp). 
After renormalization, its contributions to the effective potential (of order g 4 ) and 
to the two-derivative term (of order g 2 ) overwhelm the tree- level contributions 
when g is large* For large Yukawa coupling, therefore, quantum solitons are deter- 
mined by the fermion vacuum energy (iN/T) log det (ifl). If this has a minimum 
for given boundary conditions, the resulting configuration must have size R ~ 1/gv 
(the only scale present) and energy ~ (gv) 4 R 3 ~ d , where d is the dimension of the 
soliton. For point-like solitons (d — 0), the energy is proportional to the Yukawa 
coupling. Assuming that the large- N restriction is only a technical one to facili- 
tate calculations at strong coupling, we conjecture that these properties hold for 
quantum solitons in any strongly-coupled Yukawa theory. 
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FIGURE CAPTIONS 



1) A representative quantum kink. The solid line shows the solution of eq. (4.1) for 
the parameters y — 20 and v — 4. The dashed line shows the ansatz 
ip(z) = v tanh(z/zo), with zq given by eq. (4.6). The dotted line shows 
the classical kink, zq — 1. 

2) The function U(t). The exact fermion vacuum energy for the ansatz (p(z) = 
v tanh(^/2;o) is given by (N/x c i)yU(yz ). 
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